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Abstract. In the framework of Born-Infeld inspired gravity theories, which deviates from 
Ceneral Relativity (GR) in the high curvature regime, we discuss the viability of Cosmic 
Inflation without scalar fields. Eor energy densities higher than the new mass scale of the 
theory, a gravitating dust component is shown to generically induce an accelerated expansion 
of the Universe. Within such a simple scenario, inflation gracefully exits when the GR regime 
is recovered, but the Universe would remain matter dominated. In order to implement a 
reheating era after inflation, we then consider inflation to be driven by a mixture of unstable 
dust species decaying into radiation. Because the speed of sound gravitates within the Born- 
Infeld model under consideration, our scenario ends up being predictive on various open 
questions of the inflationary paradigm. The total number of e-folds of acceleration is given 
by the lifetime of the unstable dust components and is related to the duration of reheating. 
As a result, inflation does not last much longer than the number of e-folds of deceleration 
allowing a small spatial curvature and large scale deviations to isotropy to be observable 
today. Energy densities are self-regulated as inflation can only start for a total energy density 
less than a threshold value, again related to the species’ lifetime. Above this threshold, the 
Universe may bounce thereby avoiding a singularity. Another distinctive feature is that the 
accelerated expansion is of the super inflationary kind, namely the first Hubble flow function is 
negative. We show however that the tensor modes are never excited and the tensor-to-scalar 
ratio is always vanishing, independently of the energy scale of inflation. 
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1 Introduction 

General Relativity (GR) explains gravitational phenomena in a wide range of scales, from 
sub-millimeter to Solar System scales [1] and (if we accept the existence of dark matter and 
dark energy) even on cosmological scales from Big Bang Nucleosynthesis (BBN) time until 
today [2], Despite its observational success, attempts to modify GR have been proposed in 
the aim of alleviating some observational and theoretical problems. 

On the one hand, the dark energy problem motivates the search for modifications of the 
infrared (IR) regime of gravity in form of scalar-tensor [3-9], vector-tensor [10-18], tensor- 
tensor [19-24] and non-local [25-27] theories among others, so that the cosmic acceleration 
would be caused by a modification of GR on cosmological scales rather than induced by 
the cosmological constant or some new components [28, 29]. Also at a phenomenological 
level, some models have been put forward to account for the dark matter [30-33]. From a 
more theoretical motivation, IR modifications of gravity have been considered as potential 
mechanisms to tackle the old cosmological constant problem [34, 35], like theories exhibit¬ 
ing degravitating solutions [36-39], where the cosmological constant is effectively decoupled 
from gravity on large scales, or unimodular gravity, where the additional Weyl symmetry is 
claimed to prevent quantum corrections in the form of a cosmological constant [40-42] (see 
however [43, 44] for claims in the opposite direction). 

On the other hand, the non-renormalizability of GR also calls for modifications of grav¬ 
ity, in the ultraviolet (UV) regime this time, which should lead to a gravitational framework 
compatible with quantum physics as it is currently understood. These quantum effects are 
usually expected to regularise classical singularities present in Einstein equations. However, 
modifications of the high curvature regime have also been explored as possible mechanisms 
able to resolve geometric singularities without invoking quantum effects. From a more phe¬ 
nomenological point of view, very much like IR modifications modify the late-time cosmology 
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and, therefore, can be used to explain the current acceleration of the Universe, UV modifica¬ 
tions are expected to be at work in the early Universe. Therefore, one may wonder whether 
they could provide a new mechanism for Cosmic Inflation [45-51]. In the standard picture, 
inflation is sourced by a self-gravitating scalar field in its potential dominated regime. Among 
the currently favoured single field models [52], let us notice that some of them are modified 
gravity theories, such as the Starobinsky (or Higgs inflation) model [46, 53] which belongs to 
f{R) theories [54]. 

In this work we explore the possibility of realizing an inflationary phase through a Born- 
Infeld modihcation of gravity in the high curvature regime. The specific theory that we will 
consider here corresponds to the class of modihcations introduced in Ref. [55]. 

The original Born-Infeld theory of electrodynamics was introduced as a way of regu¬ 
larising the self-energy of point-like charged particles in electromagnetism [56]. This was 
achieved by introducing a square root structure that gives rise to an upper bound for the 
allowed electromagnetic fields. Deser and Gibbons suggested to use the same idea to resolve 
the singularities encountered in GR [57]. However, their proposal suffered from an ambiguity, 
that originated from the necessity to remove a ghost present in the metric formulation of the 
theory. In Ref. [58], the theory was considered within the Palatini formalism, putting forward 
that, in that approach, the ghost is naturally avoided. 

The phenomenological consequences and viability of Born-Infeld gravity theories have 
been extensively explored in cosmology [59-69], astrophysics [70-74], the problem of cosmic 
singularities [75, 76], black holes [77, 78], wormhole physics [79-82] and various extensions 
of the original formulation have also been considered [83-97]. Another interesting property 
of these theories is that they give specific realizations of Cardassian-like models [98, 99]. 
Recently, a natural extension of the theory was introduced in Ref. [55] where it was shown 
that the high curvature regime is free of singularity and may support a quasi-de Sitter 
expansion when the Universe is dominated by a perfect fluid having a vanishing equation 
of state (see also Ref. [91]). UV modifications of GR with additional degrees of freedom have 
been shown to provide potential candidates for dark matter [16, 100-105] so that finding 
an UV modification that may also support accelerated expansion opens the possibility of 
unifying dark matter and cosmic inflation. 

A notable topic of discussion concerning Palatini theories is the potential existence of 
anomalies and surface singularities around sharp variations of the energy and pressure densi¬ 
ties [106-111]. As argued in Ref. [112], gravitational theories containing auxiliary fields nec¬ 
essarily introduce new couplings to matter that involve derivatives of the energy-momentum 
tensor as source of the “Einstein equations”. Thus, systems in which there is a sharp change 
in the density profile may lead to strong tidal forces, although Ref. [108] has argued that 
backreaction effects might cure such a pathology. We expect these effects to appear for 
variations on the density profile of the order of the new mass scale of the theory, since the 
corrections with respect to GR are determined by such a scale. However, in Ref. [107] it is 
argued that for a specihc class of polytropic equations of state, and in a particular realization 
of Born-Infeld like theories of gravity, the appearance of the surface singularity is indepen¬ 
dent of the scale that suppresses the corrections with respect to GR. This could also be the 
case for the theory that we consider here and addressing this issue would require a detailed 
analysis that is beyond the scope of the present work. We would like to point out, however, 
that this pathology is very specihc for one theory and one class of equation of state, so such 
a problem may be avoided by assuming a high enough scale suppressing the new corrections. 
Additionally, when tidal forces start growing, additional terms in the huid description may 
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become relevant. 

In the following, we consider the Born-Infeld theory introduced in Ref. [55] and show 
that gravitating dust in the high curvature regime may be used to support Cosmic Inflation. 
The inflationary dynamics in the high curvature regime strongly depends on the speed of 
sound. In particular, getting enough e-folds of inflation to solve the usual problems of the 
Friedmann-Lemaitre model, together with having a graceful exit and a reheating era is a 
non-trivial problem. We show that the minimal model satisfying all these conditions involves 
a cascade of at least two decaying dust components ultimately decaying into a radiation 
fluid. Moreover, using Big-Bang Nucleosynthesis (BBN) constraints, such a setup gives new 
constraints on the energy scale at which the gravity modifications may take place, which ends 
up being complementary to the ones coming from astrophysical processes. The inflationary 
phase itself exhibits unique properties as for instance a bounded total number of e-folds, a 
maximal energy density and a Hubble parameter which is slowly increasing during inflation. 
Although in GR with scalar fields, such a feature would produce a blue spectral index for 
the gravitational wave spectrum, we emphasize here that tensor modes remain unamplified 
and inobservable, independently of the energy scale of inflation. 

The paper is organized as follows. In section 2 we give a brief summary of the Born- 
Infeld inspired modification of gravity that we will use as well as some of its main properties. 
Then, in section 3 we study isotropic and homogeneous cosmological solutions taking into 
account the speed of sound and show in section 4 how accelerated expansion can be obtained 
in the presence of dust. In sections 4.2 to 4.3, we show that a viable model incorporating 
a graceful exit together with a reheating before BBN requires the presence of at least two 
unstable dust components decaying one into another and ultimately into radiation. Finally, 
we discuss various attractive aspects of the model as well as some properties of the tensor 
perturbations in the conclusion. 


2 Minimal Born-Infeld theory 


The minimal extension of the Born-Infeld inspired gravity considered in the following is 
described by the action^ [55] 

5 = miMli j Tr {j 1 + mf g-^ R{T) - , (2.1) 

where Mpi is the reduced Planck mass, m\ is some energy scale (in principle unrelated to 
-^Pi); 9~^ (= 5^^) denotes the inverse of the metric tensor, R (= R^v) is the Ricci tensor 
matrix, 1 (= is the 4x4 identity matrix and Tr( ) stands for the trace operator. This 
action is treated within the Palatini formalism, i.e., the Ricci tensor is constructed out of 
an independent connection field F“j^. The factor in front of the action is chosen so that we 
recover the Einstein-Hilbert action at small curvatures and we subtract the identity inside 
the square brackets to guarantee the existence of Minkowski as vacuum solution. This can 
be seen by expanding the action at curvatures much smaller than m\'. 


S 



d^x^g^^'^R^yiV) 


1 + 0 


m{ j 


( 2 . 2 ) 


1 


Throughout this work, a hat will be used to denote matrix representation of the corresponding tensor. 
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where we see that the leading order term is nothing but Einstein-Hilbert action in the Palatini 
formalism without cosmological constant. 

Varying the action (2.1) with respect to the metric tensor yields the metric field equa¬ 
tions 

- Ti'(Af - l)mlgij,^ = (2.3) 

JWpi 

where we have defined ^ 

+ . (2.4) 

Here, we demand that 1 -|- m^‘^g~^R is a positive definite matrix on physically admissible 
solutions and we dehne M as the only positive definite matrix such that = l + m^‘^g~^R. 

The metric field equations can be written in an alternative form by using the above 
definition of the fundamental matrix in order to express the Ricci tensor as R = m\ 5 (M^ —1). 
One then obtains 

i [g{M - M-i) + (M - M-^fg 

where we have used matrix notation and the superscript “T” stands for the transposition 
operator. This equation allows, in principle, to express the matrix M in terms of the metric 
tensor and the matter content by solving an algebraic set of equations. These equations 
are, in general, non-linear and several branches may arise in the theory. Nonetheless, not 
all of them will be physical since one must require the matrix M to be positive definite. In 
addition, there is only one branch of solutions that will be continuously connected with GR 
at low curvatures (or densities). This is indeed the branch that we will choose in this work, 
although other branches can also have interesting phenomenologies [55]. 

Variations of the action (2.1) with respect to the connection T gives the remaining held 
equations in the Palatini formalism: 

Va (Vp - 2^ - SiT^pWP'^ + T^pWP'^) = 0. (2.6) 

in which we have introduced the torsion tensor and 

WP" = (M-1)'^„5"V (2.7) 

Assuming all helds to be minimally coupled to the metric, the connection equations are not 
sourced by the matter sector, i.e., the right hand side (RHS) of the connection equations 
identically vanishes^. In the following, we will be interested in torsion-free solutions and 
therefore set 7))^ = 0 from now on. Obviously one needs to check the consistency of this 
condition and we show below that it is the case for the cosmological solutions we are interested 
in. 

For vanishing torsion, taking the trace of (2.6) gives Va {y/—g = 0 such that the 

connection equations reduce to 


- Tr(M - l)g = f , (2.5) 


Va 



= 0 , 


( 2 . 8 ) 


^Some subtleties might arise when considering fermions, but we will not consider that case here. 
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for all indices. As shown in Ref. [55], if the matrix W is symmetric, an elegant way to solve 
these equations is to introduce the auxiliary metric g defined by 

= \/det (2.9) 

Plugging equations (2.7) and (2.9) into (2.8) gives 

Va = 0. (2.10) 

These equations require T to be the Levi-Civita connection associated with the auxiliary 
metric g. One can check a posteriori that having a vanishing torsion is consistent with the 
solutions where is actually symmetric. For a more detailed discussion on these points 
see Ref. [113]. 


3 Cosmological solutions 

Homogeneous and isotropic solutions associated with the action (2.1) have been discussed 
in [55] for barotropic fluids and we generalize this approach to any perfect fluid below. 

3.1 Hubble parameter 

We assume the Friedmann-Lemaitre-Robertson-Walker metric 

ds^ = —n^(t)dt^ + o(t)dx^, (3-1) 

where n{t) and a{t) are the lapse and scale factor functions respectively. The matter sector 
is modeled by a homogeneous perfect fluid so that the stress tensor and the fundamental 
matrix are assumed to be of the form 


T/ = diag[-p(t),P(t),P(t),P(t)], = diag[o(t),Mi(t),Mi(t),Mi(t)], (3.2) 


i.e., compatible with the spacetime symmetries. Since the matrix M must be positive definite, 
one has Mq > 0 and Mi > 0. The metric field equations (2.5) for the assumed background 
become 


1 

Mo 


+ 3Mi = 4 + p. 
Mo + 2Mi + ^ = 4 - P, 


(3.3) 


where we have defined the dimensionless quantities 

P 




m^M^i ’ 


P = 


P 




(3.4) 


These equations allow to obtain Mo(/5, P) and Mi(/?, P) algebraically so that we can compute 
the auxiliary metric that generates the connection in terms of the matter content. According 
to (2.9), the auxiliary metric also takes a FLRW form 


ds^ = —fP{t)dP + a^(t)dx^, 

where the auxiliary lapse and scale factor functions are given by 

a^t) = . 


(3.5) 


(3.6) 
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Since the connection is associated to the auxiliary metric g^i,, the corresponding 
time-time component of its Einstein tensor is given by the usual expression in terms of the 
auxiliary Hubble parameter H = d/a, namely 


Goo ( 5 ) = 3^2 ^ 3 


' IdMMoMi)]' 

^■4 -d^- 


(3.7) 


with H = d/a the Hubble parameter associated to the spacetime metric On the other 
hand, we can use the definition of M in (2.4) to express the Einstein tensor in terms of M 
from its definition as follows: 


Gi~g) = R{g)--lTT{r^R) 


= mlg 


- 1 - -MTr (M-M 


-1 


(3.8) 


Equations (3.7) and (3.8) will lead to the equivalent of the usual Einstein equations with a 
modified source term, since the matrix M is an algebraic function of the matter content. Eor 
the ELRW solutions, the right hand side of equation (3.8) reads 


0 * 00 ( 5 ) — 




Ml - 1 - 3MoMi + 3^ 




^0 + 9 (5 + 'f’) ^0 — 1 5 

(3.9) 


where we have used the equation (3.3) to express everything in terms of Mq only. Simi¬ 
larly, the auxiliary Hubble parameter H can be expressed in terms of Mq by noticing that 


equation (3.3) implies 


A'IqM\ 


(4 + p)Mq - 1 
3 


(3.10) 


Then, from equation (3.7) and using (3.3), one gets 


H = H 


Mq 

4[(4 + p)Mo-l] 


^ + + 


d In Mq 1 
dN \j ’ 


(3.11) 


where N = ln[a(f)] is the ”e-fold” time variable and Mq[p(N), P{N)] stands for the physical 
solution of equations (3.3). Expanding now the total derivative of Mq as 


d In Mq d In Mq dp d In Mq dP dp /'din Mq 2 ^ Mq \ 

dA^ ^ dp diV ^ dP dN ^ ^ V ^ dP ) ’ 


where = P/p is the sound speed, one obtains 


H = H 


Mq dp 

4 [(4 + p)Mo - 1] diV 


1 + (4 + p) 


/din Mq 2^1^ Mq \ 1 I 

JJr 


(3.12) 


(3.13) 


From equations (3.9) and (3.13) one finally obtains the modified Friedmann-Lemaitre equa¬ 
tion 


= 


Ml + 1{P + p)Mq - 1 


1 - 


Mq 


dp 


4[(4 + p)Mo-l] diV 


1 + (4 + p) 


, f dlnMQ „cHnMo 


2 ’ 


(3.14) 
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where H = H/m\ is the dimensionless Hubble parameter in units of the new scale mx. As 
previously advertised, our Born-Infeld inspired gravity theory within the Palatini formalism 
has led us to a version of the Friedmann-Lemaitre equation where the matter source is 
modihed. Let us stress again that Mq = Mq{p,P) so that the RHS of the above equation 
is also a function of p and P. A remarkable feature that should not be unnoticed is the 
appearance of derivatives of p and P, which will play a crucial role in the dynamics of the 
system, as we show below. Moreover, this introduces the novel effect that the background 
cosmological evolution is not only determined by the equation of state parameter w = P/p 
of the fluid as in the standard case, but also by the sound speed of the fluid Cg. 

To end this section we stress that the matter sector is assumed to be minimally coupled 
to the metric tensor g^y. In that situation, conservation of the stress-tensor yields the usual 
equation 


dp 

diV 


-3(P + p). 


(3.15) 


3.2 Physical branch 

The metric field equations (3.3) generically lead to several branches for Mq(p, P) and Mi(p, P) 
due to their non-linearity. These have been thoroughly discussed in Ref. [55] and we simply 
summarize the results here. Solving for Mq gives the cubic equation 


(4 + p)m3 + 


P(4 + p) + -(1 + p)2-4 


Mi- 


P+^i^ + p) 


Mo + - = 0. 


(3.16) 


There are three different branches of solution, but only two of them are admissible on physical 
grounds, i.e., obtained by imposing the positivity of the fundamental matrix M. Out of those 
two physical solutions, only one branch matches GR at low energy densities. Although the 
expression is not particularly illuminating, it is explicit: 


Mo = 


18(p + 4) 


2^/3 ^27V3{4 + p)V^ - A + 


2V3C 


-D 


27^3(4 + p)^/^ - A 


(3.17) 


where 


A = 3456P^ -h (108P2 -h 576P 168) - 4752^^ + (54p3 losQP^ + 1206 P 68 O) p^ 

-h (648P3 -h 3159P2 + 720P -h 2670) p^ -h (2592P^ + lOSOP^ -h 4302P -h 2616) p 
-h (72P + 96)p^ -h 9144P -h 16p® -h 1456, 

B = 144.P^ - 1248P^ -h (4P2 - 32P + 32) 2404^^ -h (9^"^ - 468.P^ -h 384.P 52) 

-h (12P3 - 80P2 - 24P -h 96) p3 -h (72P^ - 504p3 -h 256^^ - 920P -h 1024) p 

- 73456P -h 1600, 

C = (18P2 -h 144P 24) -h (144P2 -h 126P -b 200) p -b 288P^ -b (24P -b 32)p^ -b 8p^ 

- 264P -b 488, 

D = {6P + 8)p -b 24P -b 4p2 - 20. 

(3.18) 

Both Mq and M\ are well-defined, i.e., real and positive, only on a given domain in the plane 
(p, P) which has been represented in Figure 1. We have also represented in this figure various 
barotropic equations of state P = wp, with constant w, and one can single out three typical 
behaviours: 
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Figure 1. Left panel: Region in the plane (p, P) obtained by imposing that the fundamental matrix 
M is positive definite, i.e., that both Mq and Mi are real and positive valued. We see that such a 
region is basically defined by the simultaneous conditions p > — 4m|Mpj and P < m\M^y For an 
idealized barotropic fluid P = wp, with a constant equation of state parameter w, the pressure and/or 
the energy densities may end up being bounded. Right panel: We show the Hubble function as a 
function of the energy density for different equation of state parameters illustrating the 3 behaviours 
discussed in the main text. We can see that at low energy densities they all follow the usual GR 
behavior (dotted line), whereas at high energy densities {p ^ the differences appear. In 

particular, for a dust component with w = 0 we see that goes to a constant value, i.e., it gives a 
de Sitter phase. 


• For fluids with positive pressure 0 < < 1 we And that there is a maximum vaiue for 

p which is given by p < m\Mp. This is the desired property of Born-Infeid inspired 
theories, i.e., we And an upper bound for the aiiowed energy densities. 

• If 0 < tc < —2/3, there is no upper bound on p. As discussed beiow, depending on the 
behaviour of Cg, the Hubbie function can take a constant vaiue at high energy densities. 
Even though p can grow indefiniteiy, the curvature, here parametricaiiy given by H, 
remains bounded. 

• Finaiiy, for —2/3 < rc < — 1 we find that the Hubbie function grows as oc p^. In 
this case there is no realization of the Born-Infeid mechanism. Interestingiy enough, 
such a behaviour is aiso typicai of theories with extra-dimensions [114]. 

As mentioned above, these three ideal cases are only illustrative and one should keep in mind 
that in the cosmological context Cg is expected to be a function of time, and thus of the 
Hubble parameter. As a result, the trajectory followed by any realistic gravitating fluid in 
the plane (p, P) of Figure 1 is in general a curve becoming strongly non-linear as soon as p 
or P becomes of order unity. 
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On the contrary, in the low energy and pressure limit, plugging the expression for Mq 
back into the Hubble parameter (3.14) and expanding everything for small p and P, one gets 




1 _ 

+ 


- 1 5 

——Pp + 





(3.19) 


As expected, this expression matches the usual Friedmann-Lemaitre equation for p <C m|Mpj 
and P <C As a result, the energy scale at which one should expect deviations from 

the standard GR case is determined by the geometrical mean of m\ and Mpi. 


4 Inflationary scenario 

In the previous section we have reviewed the cosmological evolution for the Born-Infeld 
inspired gravity theory under consideration. We have shown that for a perfect fluid with 
equation of state parameter satisfying —2/3 < w < 0, the energy density may not be bounded 
from above, but the curvature is. In particular, as we show in more detail in the next 
section, for a dust gravitating fluid having w = 0, the Hubble parameter becomes nearly 
constant at high energy densities thereby allowing a quasi-de Sitter expansion typical of 
an inflationary epoch [115]. Since the theory matches GR at low energies, the inflationary 
graceful exit is always naturally realised within our Born-Infeld inspired gravity when such 
a regime is reached. However, with only dust, the Universe would end up being matter 
dominated after inflation and, in order to produce a radiation dominated Universe, one needs 
to implement some mechanism allowing the dust to decay into radiation. This is analogous 
to the reheating period within the standard inflationary picture where the inflaton decays at 
the end of inflation giving rise to radiation made out of relativistic degrees of freedom. 

Thus, in section 4.2, instead of stable dust, we will consider an unstable dust fluid 
decaying into radiation. However, in this scenario the sound speed is no longer vanishing, 
see equation (3.14), and this has a crucial effect since now the existence of a quasi-de Sitter 
era becomes possible only within a finite range of energy densities. Although inflation can be 
made to last long enough by adequately fixing the decay rate of the dust component, we show 
that the duration of the reheating era ends up being necessarily longer than the inflationary 
period. On the other hand, the decay rate will also determine the beginning of the radiation 
era, so the duration of inflation and the reheating period are linked. This in fact prevents 
the model from solving the flatness problem of the standard FLRW cosmology. 

In section 4.3, this issue is solved by considering a cascade of decaying dust fluids which 
ultimately end into radiation. In such a case, inflation can be realized with a graceful exit 
onto a reheating era. The solution comes about because while the duration of inflation 
depends on the whole cascade of dust fluids, the reheating duration depends only on the 
decay rate of the last component thereby making the model viable. 


4.1 Stable dust 


Let us start by considering the gravitating fluid to be pure dust and conserved 

by 


p 

w = — = 0, 
P 


Cs = 0, 


dp 

dA^ 


= -3p. 


characterized 

(4.1) 


Plugging these conditions into equations (3.14) and (3.17) gives the Hubble parameter as 
an algebraic function of p only. If we take the high energy density limit p ^ 1 (i.e., in the 


- 9 - 





Figure 2. Hubble parameter square as a function of the stable dust energy density p (left panel). 
In the Born-Infeld regime, p 3> 1, the Hubble parameter is nearly constant and slowly grows H > 0 
to generate a superinflationary era. The red dotted curve corresponds to the asymptotic expansion 
of equation (4.2). In the right panel, the first Hubble flow function ei is represented as a function of 
p. In the superinflationary regime ei < 0. 


Born-Infeld regime), we find 


H = 


^ 8 ^ 3 - 3^6 ^ ^ 




(4.2) 


so that the Hubble parameter becomes nearly constant and the expansion of the Universe is 
accelerated. The exact dependence of H with respect to p as well as the above expansion 
are represented in Figure 2. Strictly speaking, acceleration of the scale factor, i.e. inflation, 
occurs as long as the first Hubble flow function ei = —dln(H)/dA^ is less than unity. From 
equation (4.2) we obtain 


ei = -- 


2p 


+ 0 


such that for p !:$> 1 one has ei < 0. Negative values of ei cannot be obtained within General 
Relativity and correspond to superinflation [116]. Given that the gravitational sector is 
not described by the usual Einstein-Hilbert term, let us emphasize that it is not possible to 
deduce from ei < 0 that the primordial perturbations will have a nearly scale invariant power 
spectrum with a blue spectral index. The cosmological perturbations are discussed in more 
detail in the conclusion where we show that tensor perturbations are in fact not generated. 

Finally, as can be seen in Figure 2, as soon as p < 1, the Hubble parameter evolution 
becomes nearly identical to GR (up to some relaxation oscillations) and inflation naturally 
ends. However, in the subsequent decelerated expansion, the Universe is and will remain 
matter dominated, which is in contradiction with both BBN and the existence of the Gosmic 
Microwave Background (CMB). The most natural extension to produce a radiation domi¬ 
nated Universe after inflation is to assume that the dust component is actually decaying into 
radiation. We explore this possibility in the next section. 


4.2 Decaying dust and radiation 

In order to generate a radiation dominated universe after inflation, let us now consider an 
unstable dust component that decays into radiation at a constant rate F. Thus, the matter 
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sector consists of two fluids, decaying dust and radiation, in interaction, which are described 
by the following coupled equations: 

^ + 3///ii = -r/ii, ^ + 4/7p, = rpi, (4.4) 

where pi and pr are the energy densities of dust and radiation respectively and H is given by 
equation (3.14). Again, since the matter sector is minimally coupled, we have conservation 
of the total energy density p = p^ + pi. However, the total pressure does no longer vanish and 
reads P = Pr/3. In addition, we are also in the presence of a time-dependent sound speed Cg. 

We can rewrite the equations of the interacting dust and radiation system in a more 
convenient manner by introducing the radiation fraction 


X = El = 

~ p pr +pi' 

so that equation (4.4) can be recast into equations of evolution for p and X 
^ = -(3 + X),. ^ = + 

where F = F /mx and P = Xpl‘P. From the definition of Cg = P/p, one gets 

2 _ 4A 1 - X f 

“ 3(3+ X) “ 3(3 + X)H ■ 


(4.5) 


(4.6) 


(4.7) 


Because the above expression explicitly involves H, and equation (3.14) depends on Cg, the 
resulting evolution of the system is highly non-linear and crucially depends on the functional 
form of Cg. As discussed before, because the speed of sound “gravitates” within Palatini the¬ 
ories, this is expected and emphasizes the importance of considering more than the equation 
of state parameter in the fluid description in order to obtain the background cosmological 
evolution. 

Plugging equation (4.7) into equation (3.14) yields an algebraic equation for H that can 
be analytically solved in terms of P and p as^ 


H = 




1 + 


HP + p) 


Mo 


(4 + p)Mo - 1 


1 + (4 + p) 


(4-hp)cs^F 9 In Mo 
(4 + p)Mo-l dP 
f (9In Mo ^ ^dlnMo'^ ’ 

V 9/0 ^ dP ) 


(4.8) 


where we have defined 

2_ 4X 2_ 1-^ 

3(3+ X)’ 3(3+ X)' 


(4.9) 


This expression matches equation (3.14) if one sets both the decay rate F and the radiation 
fraction X to zero, as one may expect. However, for non-vanishing F, the limit X —)■ 0 does 
not give back the pure dust behaviour as there will always be a term proportional to F in 
the expression of H. This is important for the initial conditions of the inflationary scenario 


^When solving this equation we obtain two branches. We choose the one corresponding to an expanding 
rather than contracting universe. 


- 11 - 
















since it will set a maximum value for the total energy density. More precisely, expanding 
equation (4.8) at large p and small X we obtain for the initial value of the Hubble parameter 


= \ - - 


21 + V2- r 


24 


r - 


6^2 


P + 


3^3 




288 




/ -124 + 189^2 -8548 + 4217^2 ^ 


12^3 


1152 


X + 8\/3 - 3^6 + 


72 

248 - 249^21 
64 i 


+ 


- 59 V6 + rU + 0 f 1, 

3^3 4608 j p \p^' 


(4.10) 

As opposed to pure stable dust, there is now a maximal value of p at which the Hubble 
parameter vanishes. Such a behaviour happens for any value of X, even for A = 0, and 
therefore is not due to non-vanishing pressure but induced by a non-vanishing speed of 
sound. Let us stress that this is only relevant for the initial conditions since, if we have a 
non-vanishing value of L, the evolution of the system will immediately generate a certain 
amount of radiation. Under these considerations and at leading order in small X and L, one 
gets H (pmax) = 0 for _ ^ 

/ f 4 


Pmax — 


=x 


3 V 6 V 2 3^3 

where we see that for X = 0, we obtain the non-vanishing maximum value 


Pma,x — 


24 

7^ 


(4.11) 


(4.12) 


This condition guarantees that the universe starts in an expanding phase, so the inflationary 
period will be achieved. We should emphasize that this bound appears if we impose a 
vanishing radiation component initially. If we have some radiation initially, then (4.11) 
should be used instead. Interestingly, initial energy densities higher than pmax continuously 
map into negative values of the initial Hubble parameter. This suggests that some of our 
inflationary solutions could be extended back in time with a bounce [117-119]. 

In order to discuss the inflationary phase only, we will assume in the following that the 
Universe starts its evolution with p(0) = pini < Pmax, in the regime in which H is independent 
of p, and with only decaying dust, i.e. having A(0) = 0. Equations (4.6) can be approximated 

by 


dX 

dX 


-A- + I. 

H ’ 



(4.13) 


the solutions of which are 



p{N) = pi„ie 


(4.14) 


Therefore, while the Universe is inflating, the radiation fraction remains constant and given 
by T/H. This in turn allows to obtain the condition guaranteeing that the system remains 
inside the physical region of the theory as depicted in Figure 1. The radiation component 
gives rise to a positive pressure given by P ~ ^Xp so that, according to the above solution 
for X, we have P ~ Hp/{3T) so the condition to be inside the physical region P < 1 leads to 
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a bound for p given by p < \/8/(3\/3f). It is of the same order of magnitude as the bound 
given by (4.12), although obtained from a different condition. 

On the other hand, the total energy density is driven by the dust component and 
decreases exponentially fast. Inflation ends for ei{p) = 1, i.e. when we exit the Born-Infeld 
regime and the evolution equations become close to those in GR. In order to obtain the 
duration of the inflationary era, we only need an order of magnitude estimate for p(A^end) 
which can be obtained by solving H{p) = p/3 by using the asymptotic expansion (4.10). 

At leading order, one gets p(Aend) — 8 and, therefore, the maximum number of e-folds of 
inflation is 




(4.15) 


In order to have more than 60 e-folds of inflation (which is the typically required minimum 
duration of inflation), the decay rate of the dust component should be very small, namely 
r < 10“^®mA, which compromises the appearance of the radiation era after inflation. For m\ 
of the order of the Planck mass, this implies that the dust component should have a lifetime 
around 10^^ s, i.e. much longer than the (GR) age of the Universe ~ 10^^ s. 

Considering a mass scale nix much larger than the Planck mass does not help either. 
Indeed, the reheating proceeds only when the cosmological evolution matches GR, and the 
radiation era starts when the reheating is completed at X(Areh) — 1- Assuming X <C 1 
during reheating, we can approximate the evolution equations by 


^ PiN > iVend) - P{Nend) 6 ^ 8e , (4.16) 

where AN = N — A^end- Therefore, with A(Aiend) — 0, one gets 

> iV..,) ^ r (4,17) 

pT^n) 5V2 

and solving for X ~ 1 gives 

AiV,eh ^ - In = 2AiVinf + -ln(^-y (4.18) 


One concludes that, independently of f and mx, the reheating era typically lasts twice the 
number of e-folds of inflation. Because reheating is a decelerating era, this catastrophically 
prevents the inflationary era to solve the flatness and horizon problem. The problem is that 
the decay rate determines both the duration of inflation and the reheating period so that 
both are intimately linked. In the next section we avoid this difficulty by introducing more 
than one decaying dust component. 


4.3 Cascading dust and radiation 

In order to achieve a viable inflationary scenario leading to a radiation dominated universe, 
one is led to consider a cascade of unstable dust components decaying one into another and 
ultimately producing radiation. Intuitively, the number of e-folds required to reheat the 
Universe will be set by the lifetime of the most stable specie whereas the speed of sound 
evolution, and thus inflation, is expected to be sensitive to all species or, equivalently, all the 
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unstable components will contribute to generating a radiation fraction and, thus, to the total 
pressure. For this scenario, one finds that the relation between AA^reh and AAjnf is indeed 
relaxed. The conservation equations for each species are given by 


^ + 2,Hpi = Ti_ipi_i - TiPi 

at 

% + m-p, = TnPn, 


i = 1,..., n 


(4.19) 


with Fq = po = 0 and where pi are the energy densities of the dust components, Fj are 
the decay rates of the corresponding particles and p^ is the energy density of the radiation 
component as before. Introducing the relative fractions 




Pi 


E n - I — 

i=l Pi + Pv 


X = 


Pr 


E n - I - 

i=l Pi + Pr 


one has now to solve the coupled system of equations 


(4.20) 


dX, 

diV 



Xi + 


F.-i 

H 


Xi-i 


dX 

dX 


(X - 1)X + ^X„, 


(4.21) 


with H given in equation (3.14). Again, the total energy density p is conserved so that 
it satisfies (4.6). The total pressnre is still driven by the radiation component and reads 
P = Xpj3 whereas the sound speed becomes 


c 


2 

S 


4X Xn fn 

3(3+ X) 3(3 + X)Fr' 


(4.22) 


4.3.1 Inflationary regime 

Comparing this expression to equation (4.7), one deduces that the Hubble parameter for 
cascading dust and radiation is given by equation (4.8) with 


2 _ 4X 2 _ 

3(3+ X)’ 3(3+ X) 

As a result, the Hubble parameter has an explicit dependence only on the total energy density 
p, the radiation fraction X, the fraction X„ of the last specie in the decaying cascade as well 
as its associated decay rate F^ (into radiation). Expanding the Hubble parameter at large 
value of p and small values of X and X„, one gets 


(4.23) 


u — / “ 21 + \/2 ^ ^ r„x„ _ ^ 

'^3 24 6 ^ 2 ^ 


4 + 173^2 


3^3 


288 


F„X„ Xp - 


4 +V2f, 


72 


FnX^Xp 


-2 


+ 


+ 


-124 + 189^2 -9556 + 4169v^- 


12^3 


+ 


1152 


F^X^ I X + I 8\/3 — 3\/6 + 


248 - 249^2 


64 


^3 (856 + 493v^) 7 (67488 - 44323v^) - 


90 + 63^2 


4608 


FnX^ 


^ + o( ^,x^x, 


FnX, I - 


(4.24) 
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Up to some numerical coefficients, this expression is formally identical to equation (4.8) 
with the replacement F —)• Therefore, starting within the superinflationary regime at 

vanishing radiation now requires p < Pmax with 


24 


Pmax 




(4.25) 


and the duration of inflation can, a priori, be made longer by having small values of 
instead of F^. However, Xn is not arbitrary but given by solving the set of equations (4.21). 
The crucial point to realize here is that, in the previous case with only one dust component, 
having small X required having small F as well, i.e., the two conditions were related, whereas 
in the cascading case the condition will depend on the whole set of decay rates, as we show 
in the following. 

Assuming H to be almost constant with p(0) = < Pmax together with 2fi(0) = 1, 

Aj^i(O) = 0 and X(0) = 0, one can find an approximate solution for all the species. One 


gets 


dXi 

d/V 


-—All, 

H ’ 


dXi 

diV 



F.-i 

H 


Xi-i. 


(4.26) 


The first dust component can be immediately integrated as 


Ai(Af) = exp () ~ 1, 


(4.27) 


where the last approximation requires N F7/Fi. For instance, if we require the duration 
of inflation to be O(lO^), we need Fi < 10^77. Obviously this is nothing but the condition 
that the dust component should remain stable throughout the whole inflationary phase. The 
remaining components of the cascade Xi’s are given by the quadrature 




■p. rN 

(N) ~ f 
Jo 


Plugging Xi ~ 1 into this equation for i = 2 gives 


X2{N) ~ 

F 2 L i H 


(4.28) 


(4.29) 


where the last approximation is valid for N <C .H/max(fi,f 2 ) (which is again related to the 
stability condition of the cascade during inflation). From this expression for X 2 {N) one can 
compute X 3 (X) and so on such that the full hierarchy is given by 


XniN) ~ 


?\rn —1 ) 

(^n 

i=i 


tj_ 

H 


(4.30) 


again for N <C 
cascade of dust 
given by 


77/max({Fj}j=i^„_i). As a result, concerning the inflationary phase, the 
components behave as one unstable dust fluid with an effective decay rate 


F = FnXn 


nr=i r. 

Hn-l 


~ (n-1)! 



(4.31) 
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Although f is no longer constant, its dependence on the e-fold number N is only a power 
law such that pmax acquires the same dependence on N. Since the total energy density p{N) 
decreases much faster, as exp(—3A^), it is only necessary to require p < pmax at the beginning 
of inflation. In order to show this more precisely, we need to impose /o/pmax ^ 1 throughout 
the whole inflationary period. Then, we only need to compute the number of e-fold at which 
p/Pmax ^ takes its maximum value'^. We find that this happens for N = {n — l)/3. Since 
we expect inflation to last for about 60 e-folds, the obtained value will only make sense if it 
is smaller than ~ 60. This is indeed the case if n — 1 < 180, so we will assume this in the 
following. Thus, taking N = (n — l)/3 m the previous expression gives the maximal number 
of inflationary e-folding 


AWnf - 3 111 


V3 


nr=i r* 


n — 1 , 

H-In 


24 1 , 


(4.32) 


For instance, taking all the fj equals, getting 60 e-folds of inflation reqnires Tj < 

for a mixture of n = 3 dust components. Taking mx = 0{Mpi), the minimum lifetime of 

each specie is required to be greater than 10“^® s. 


4.3.2 Reheating 

The reheating can be dealt as in the previous section. Provided the radiation component re¬ 
mains snbdominant, X <C 1, one has H ~ x/ p/3 with p{N) ~ p(A'end) exp(—3AAi). Plugging 
these expressions into equation (4.26) gives the evolution of the density fraction of all the 
species during reheating, up to the end at which the radiation fraction becomes significant. 
For Xi, one gets 


Xi{N > Nend) = Ai(iVend) exp 


_2^J_ / 3AN/2 _ ^ 

3 Rend ^ 


(4.33) 


where Rend = -\/p(Rend)/3 — x/3/3. Moreover, under the reasonable assumption that Xi 
does not evolve much during inflation, i.e. AWnf <C f y^ 8/3, one has Ri(Rend) — 1- 
During reheating, Xi[N) remains almost constant as long as AR <C ARi where 


AW = W - Aend = ^In ^ 


(4.34) 


For the other species, one obtains for i > 1 

2 f. 


W(R)~W (Rend) exp 




liz 

Hend 


f-N 

- / exp 
d -lA^end 


3 Rend 
2 Vi 


3AN/2 _ 

(^3Ay/2 _ 3AN/2\ 3^ 

3 Rend V J 2 ^ 


(4.35) 


Xi_i{y)dy. 


From the inflationary solution, Rj(Rend) is given by equation (4.30). The first line in the 
previous expression encodes the evolution in the absence of sources, such that for AN < 

Xi keeps its value acquired during inflation. For AN > ARj, the source-free evolution is 
decaying exponentially of exponentially fast snch that Xi is completely driven by the sonrce 

"^As in the previous cases, this is quantitatively similar to imposing the total pressure to be smaller than 
i.e., P = l/3Xp = l/3f„A„p < 1. 
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given in the second line of equation (4.35). As a result, either Xi takes its value at the end 
of inflation, or the source term drives its dynamics. For this reason, we only consider the 
source term from now on. 

For X 2 , plugging equation (4.33) into (4.35), one has 


A 2 (iV) ~ ^ f"" exp (I ^ 3 _ (4 gg) 

-Hend J N^nd V-J -^end ^ / 

The argument under the integral emphasizes the typical time scale AA^i 2 involved in the 
evolution of X 2 sourced by Xi. It is defined by the transcendental equation 



2|r2-fi 
3 Send 


„3AAfi2/2 
^ ? 


(4.37) 


whose solution can be expressed in terms of the Lambert function 


ANu = 



(AlilziA) 3 _^end \ ^ 2 ^ \ 

3 77end ) 3 \^2|r2-ri|y 3 \^|r2-ri|y’ 


(4.38) 


the asymptotic limit holding for |r 2 — Fil ^ Lfend- For 1 < AN < AN 12 , we obtain the 
approximate solution 


^ 2(1 < AN < AN 12 ) ~ 


2 Fi 


3N 


exp 


end 


/ 2 r2 ^3AAf/2 
V 3 .^end 


+ 5AiV 


(4.39) 


showing that X 2 first exponentially grows for AN < AN 2 and then disappears as a double 
decaying exponential for AN > AN 2 . However, such an evolution may be interrupted before 
completion if AA 12 < AN 2 . 

For AN > AA 12 , one has to distinguish whether r 2 <riorr 2 >ri. In the former 
situation, the argument of the exponential under the integral of equation (4.36) becomes 
strongly negative and the integral saturates for Ay ~ AA 12 . One has 


X 2 {AN > ANu) ^ exp + ^AAi 2 ) , r 2 < f 1 . (4.40) 

This expression shows that for the “interrupted” case, AA 12 < AN 2 , the exponential growth 
of X 2 is interrupted at AiVi 2 and X 2 remains constant up to AA 2 at which it hnally disap¬ 
pears. From equation (4.38), the stationary (and maximal) value of X 2 is 


X 2 {ANi 2 < an < AN 2 ) ~ 


Fi 

f 1 - r2 ■ 


(4.41) 


Under the same hypothesis, AN > AA' 12 , let us discuss the case f 2 > fi. The argument of 
the exponential in equation (4.36) becomes positive and the integral blows up as a double 
exponential such that one can neglect the term in 3Ay/2. The integral is now given by an 
exponential integral function such that 

JAfgnd V-J -"end / 6 H^nd J 
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Figure 3. Evolution of the relative dust fraction Xi (green) and radiation X (red) for three dust 
components during reheating. One may notice the transient stationary behaviour described by equa¬ 
tion (4.41). Reheating ends at the time the last dust component of the cascade decays. 


Taking the large argument limit in the previous expression one hnally gets 

X2(AiV> AiVi2)~:^^^expf-^-^e3^'^/2A f2 > fi. (4.43) 
t 2 — t 1 V -Hend / 

As one may expect, for r 2 >ri, the late time evolution of X 2 becomes completely driven by 
the evolution of Xi and thus X 2 disappears when A A ~ AAi. 

Equations (4.39), (4.40) and (4.43) give a good approximation of the evolution of X 2 
at all stages. Plugging back these expressions into equation (4.35) would give in the same 
way the evolution of all Xi by recurrence. In fact, we do not need to perform any additional 
calculations. Ultimately, either r 2 < Pi and X 2 decays as in equation (4.40), i.e. completely 
disappears for A A > AA 2 ; or r 2 > fi and A 2 decays as in equation (4.43), i.e. for A A > 
AAi. Both equations (4.40) and (4.43) are of the same functional form as equation (4.33) 
such that the ultimate behaviour of all the Xi will be of the same functional form. As 
a result, after some transient evolution, the component Xi is expected to disappear when 
A A ~ sup{AAj,j < i}. 

In figure 3, we have numerically solved the evolution equations for n = 3 decaying dust 
components in addition to radiation having Tj+i < Tj. As can be seen in this plot, the 
expected behaviours of Xi are recovered. While Xi decays, Aj>i exponentially grows to 
reach a stationary regime and finally disappears. Reheating ends when the Universe contains 
a significant amount of radiation, namely for A(Ai.eh) — 1- Because 


X = (4.44) 

i=l 


the end of reheating corresponds to the time at which the last dust component disappears, 

i.e. 


AAreh 


max (AAj) ~ - In 
3 


V6 

min(rj) 


(4.45) 
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Figure 4. In the left panel we show the lower bound on mx from the constraint given in (4.49) as a 
function of the number of dust components n when they share the same decay rate Fi. The allowed 
region is above the curve so we see that only for n = 1 the constraint can not be made below the 
Planck scale. The lower bound (4.50) is represented by the dotted line and we see that is only satisfied 
for n < 20. In the right panel we show the different bounds in the plane (T^jm^). The orange region 
is the bound on F^ for reheating to end before BBN. The blue region is excluded from the fact that the 
dust components need to be stable during inflation as expressed in (4.50). Finally, the green curves 
correspond to bounds on mx in order to have, at least, AiVinf — 60 e-folds of inflation. 


Notice that only the most stable component contributes to the duration of the reheating 
period, whereas, as can be seen in equation (4.32), the duration of inflation depends on the 
whole cascade. This is the crucial feature that allows to make the model viable by introducing 
the cascade. One can readily understand that with only one dust component, both durations 
are linked and this was precisely the obstruction to make the model viable with only one 
dust component. 

4.3.3 BBN constraint 

Comparing equations (4.32) and (4.45), one can obtain a long enough inflationary era to 
solve the usual Big-Bang problem without having a too long reheating era. Reheating should 
however end before the onset of Big-Bang Nucleosynthesis, i.e. for preh > Pnuc- According to 
the previous section, 

Preh = ~ ^ [min(f j)] ^ , (4.46) 

such that the BBN bound for pnuc — 10 MeV reads 


min(R) > ^ 10-^1 Mpi. (4.47) 

For three fluids having the same decay rate, and mx = O(Mpi), getting 60 e-folds of inflation 
were requiring only Fj < 10“^® mx so that the BBN constraint is trivially satisfied for that 
case. 


- 19 - 















In the case where all the Fj are equal, one can use the BBN bound to actually constrain 
the acceptable values of mx in order to have a viable inflationary model. From equation (4.32) , 
one has 

SAA^inf ~ —nln (^ In 
\mxj 2 

If we saturate the minimal value for F* we can obtain the following lower bound for mx by 
requiring that inflation should last longer than ~ 60 e-folds: 


24 


(n — 1)^ 


1 n 

+ ln[(n-l)!] + ^. (4.48) 


mx > 


n — 1 

r ^ 1 

1 /n 

\/3/^nuc n 1 

r V8 1 

Vu 

(n — l)(n — 1)! 

2Mpi 

(n — l)(n — 1)! 


-94+3A]Vi„f 

e ^ 


(4.49) 

This bound implies that for one single dust component n = 1 one would need mx > lO^^Mpi, 
well above the Planck mass. In addition to be unnatural, we have shown earlier that the 
duration of reheating would be twice longer than inflation in this case. Already for n = 2 we 
however obtain the bound mx > 10“^Mpi, which is below the Planck scale. Finally, in the 


limit of many intermediate species (but smaller than 180 as we required before) the smallest 
possible bound is mx > 10“^^Mpi. There is an additional constraint that we have on mx 
coming from the stability of the cascade during inflation, i.e., the decay rates should be smaller 
than the expansion rate during inflation for the dust components to remain throughout the 
inflationary phase. This was already stated as the condition AAjnf H/Ti in order to have 
Xi ~ 1. Thus, if we again saturate the value of Fj, we obtain the additional bound 


mv > y|Fi AiVinf ~ 10-39Mpi, (4.50) 

where we have used AAjnf ~ 60. We see that this bound is more stringent than the one 
obtained in the large n limit. This means, that there will be a maximum value for the 
allowed number of species such that (4.50) is satisfied. Such a maximum value is n < 20 as 
can be clearly seen in Figure 4. In that figure we have also reported all the bounds in the 
plane {Ti,mx). 


5 Conclusion 

In this paper, we have discussed a scenario in which inflation is a natural outcome of a 
modification of gravity in the high curvature regime. 

By considering minimal Born-Infeld gravity in the metric-affine formalism, we have 
shown that accelerated expansion of the Universe can be induced by gravitating dust at energy 
densities higher than mx being the new mass scale of the theory. At lower curvatures, 

or energy-densities, GR is recovered such that these theories automatically implement a 
graceful exit of the accelerated phase without the need of scalar fields (and/or negative 
pressure fluids). Graceful exit is however not enough to ensure a viable inflationary scenario 
as the Universe should reheat to become radiation dominated after inflation, and such a 
reheating era should end before BBN. 

Reheating within our scenario can be implemented by considering that inflation in the 
Born-Infeld regime is supported by a chain of unstable dust components that ultimately 
decay into radiation when the GR regime is recovered. Because in the Palatini formalism 
“the speed of sound gravitates” or, in other words, the background evolution not only depends 
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on the equation of state parameter, but also on the sound speed, we have shown that such a 
scenario is severely constrained thereby making the model relatively predictive. For instance, 
the total number of e-folds cannot be made arbitrarily large, see equation (4.32), and is 
completely fixed by the decay rate of the dust species. As a result, observing a tiny but non¬ 
vanishing spatial curvature today would be a natural outcome of such a scenario, unlike in 
most inflationary scenarios where the typical duration is usually much larger than O(lO^) e- 
folds and, therefore, any initial spatial curvature is heavily suppressed. Another requirement 
of a successful reheating phase is that at least two decaying dust components should coexist 
during inflation, cascading one into another. Indeed, when only one unstable dust specie 
decays into radiation, both the duration of reheating and inflation are completely fixed by 
the dust lifetime. In that situation, one obtains that reheating lasts twice the duration 
of inflation, which is in tension with cosmological observations. For more than two dust 
components, the model works provided the new mass scale mx is not too low (see figure 4), 
although the constraints are fairly mild and having m\ > 10“^®Mpi would suffice. 

Interestingly, due to the existence of a non-vanishing speed of sound Cg during inflation, 
there is a maximum energy density Pmax) given by equation (4.12), above which the infla¬ 
tionary phase cannot exist anymore and could be superseded by a bounce. In addition, the 
presence of a positive pressure given by the radiation component also sets an independent 
upper bound on the energy density in our Born-Infeld inspired modification of gravity. As a 
result, the energy densities of the (unstable) dust components are self-regulated and cannot 
take arbitrarily large values in the past. On the one hand, cascading dust inflation provides 
a framework in which, at the background level, various criticisms of the standard inflationary 
scenario, and criticisms of the alternative models, seem to be addressed. On the other hand, 
there are still various points to be clarified. For instance, the nature of the dust components 
has not being specified and could range from super-heavy or super-cold dark matter particles 
to black holes. Being more specific on the origin of the gravitating dust should put additional 
constraints on the acceptable decaying rates, and thus on the duration of inflation and reheat¬ 
ing. It will be also interesting to discuss if these particles can play the role of dark matter. 
Most importantly, we have not derived the primordial power spectra that would source the 
cosmological perturbations. This is indeed a non-trivial problem within the Born-Infeld class 
of theories as we illustrate now. 

A distinctive feature of the model is that it exhibits superinflation, namely a negative 
first Hubble flow function ei < 0. Within GR, one would immediately conclude that the 
spectral index of the gravitational waves is blue, and from equation (4.3), that the second 
Hubble flow function €2^3 such that the spectral index of the scalar modes would be super- 
red. However, this is not correct as in the Born-Infeld regime the perturbations evolve in a 
completely different manner than in GR. The problem for the tensor modes hij = a~^5gij 
has been addressed in Ref. [120] for a large class of Palatini theories and it was shown that, 
in the absence of anisotropic stresses (as expected in our inflationary scenario), they verify 


hij — hij , 


-// I 

hij + 




hij — 0, 


(5.1) 


where hij = a~‘^5gij stands for the tensor perturbations of the auxiliary metric and jlij = ahij 
is the usual quantized mode function. In the above equation a “prime” denotes derivative 
with respect to the auxiliary conformal time f) (defined by h = a). The mode evolution for 
Jlij is identical to the one of GR, but with respect to the auxiliary metric and coordinates. 
For a quasi-de Sitter expansion H ~ in the metric g, one gets a{fj) oc i.e. 


- 21 - 


kination in the auxiliary metric. As a result, the effective mass in (5.1) reads —a"la = 
+ \/{A'rf‘) > 0 and no amplification of the tensor modes occurs. As a result, and in spite of 
the superinflationary nature of the cascading dust model, no primordial gravitational waves 
are generated, independently of the energy scales involved. This is indeed a very distinctive 
feature of this inflationary model. In particular, any detection of S-modes in the CMB 
generated by primordial gravitational waves would rule out the model. 

The situation is more complex for the scalar modes and, at the time of this writing, no 
conclusion can be drawn on their viability. Their equations of motion will now involve the 
scalar sources, namely perturbations in the dust components and one would need to quantize 
the fluid degrees of freedom directly. This could be done either by resorting to an effective 
theory for fluids [121-125] or by computing the spectrum of thermal fluctuations. We leave 
however the explicit computation for a future work. 
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